Universal shape law of stochastic supercritical bifurcations: Theory and experiments. 
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A universal law for the supercritical bifurcation shape of transverse one-dimensional (ID) systems 
in presence of additive noise is given. The stochastic Langevin equation of such systems is solved 
by using a Fokker-Planck equation leading to the expression for the most probable amplitude of 
the critical mode. From this universal expression, the shape of the bifurcation, its location and 
its evolution with the noise level are completely denned. Experimental results obtained for a ID 
transverse Kerr-like slice subjected to optical feedback are in excellent agreement. 
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In nature, most of physical systems are subjected to 
fluctuations. For a long time, the effects of these fluctu- 
ations were either considered as a nuisance (degradation 
of the signal-to-noise ratio) or ignored because it was 
not known how to handle them. Since three decades, 
a wealth of theoretical and experimental researches have 
shown that fluctuations can have rather surprisingly con- 
structive and counterintuitive effects in many physical 
systems and that they can be figured out with the help 
of different analysis tools. These situations occur when 
there are mechanisms of noise amplification or when noise 
interacts with nonlinearities or driving forces on the sys- 
tem. The most well-know examples in zero dimensional 
systems are noise induced transition []J and stochastic 
resonance [2j. More recently, examples in spatially ex- 
tended system are noise induced phase transition, noise- 
induced patterns (see [3j and references therein), noise- 
sustained structures in convective instability stochas- 
tic spatio-temporal intermittency 0], noise-induced trav- 
elling waves jo], noise-induced ordering transition 0] and 
front propagation A direct consequence of noise ef- 
fects is the modification of the deterministic bifurcation 
shapes. The critical points and the physical mechanisms 
are masked by fluctuations. It is important to remark 
that the critical points generically represent a change of 
balance between forces. Hence, the characterization of 
noisy bifurcations is a fundamental problem due to the 
ubiquitous nature of bifurcations. For instance, the su- 
percritical bifurcations transform into smooth transitions 
between the two states and the subcritical bifurcations 
experience hysteresis size modifications. In absence of 
noise, the shape of a bifurcation and its characteristics 
are given by the analytical solution of the deterministic 
amplitude equation of the critical mode jjj. On the other 
hand, in presence of noise, no such analytical expression 
can be obtained from the stochastic amplitude equation. 
In this latter situation, the below and above bifurcation 
point regimes can be treated separately, but without con- 
tinuity between their respective solutions. For instance, 
in noisy spatially extended systems in which the systems 



are characterized by the appearance of pattern precursors 
below the bifurcation point and by established patterns 
that fluctuate above this point [lfj, the precursor ampli- 
tude [111], obtained from the linear study of the stochastic 
equation, diverges at the bifurcation point and does not 
connect to the "mean" amplitude of the fluctuating pat- 
tern, obtained from the deterministic equation. To our 
knowledge, no universal analytical expression of the crit- 
ical mode amplitude, describing the complete transition 
from below to above the bifurcation point, exists for the 
supercritical bifurcations in presence of noise. 

In this letter, we propose a universal description of the 
supercritical bifurcation shapes of ID transverse systems 
(either uniforms or very slowly varying in space) in pres- 
ence of noise that is also valid for the second order bi- 
furcations of temporal (zero-dimensional) systems. More 
precisely, we give a unified expression for the most prob- 
able amplitude describing the supercritical bifurcations in 
presence of noise, including the noise level and the bi- 
furcation point location. The systems under study are 
described by stochastic partial differential equations of 
the Langevin type [T2| (first order in time and with lin- 
ear noise terms) involving additive white noise. Firstly, 
we solve the Langevin equation describing the stochas- 
tic dynamics by using a Fokker-Planck equation for the 
probability density of the critical mode amplitude. Then, 
from the stationary distribution of this amplitude, we de- 
duce the bifurcation shape by means of the most probable 
value of the pattern amplitude. Finally, the comparison 
with experimental results obtained in a Kerr-like slice 
subjected to ID optical feedback is given and leads to an 
excellent agreement. 

Let us consider a ID extended system that exhibits a 
supercritical spatial bifurcation described by 



d t u = f(u,d x ,{n}) + y/rjCfat) 



(1) 



where u (x, t) is a field that describes the system under 
study, / is the vector field, {//} is a set of parameters that 
characterizes the system, r\ is the noise level intensity, and 
£(x, t) is a white Gaussian noise with zero mean value and 
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Figure 1: Averaged Fourier transform modulus S(k) of the 
field u for the supercritical Swift-Hohenberg model [3] ob- 
tained below the bifurcation point. The inset figure corre- 
sponds to an instantaneous snapshot of S(k,t). 

correlation (( l {x, t)( j (x, €)) = 6 ij S (f- t) 5 (x- x). 

We assume that the associated deterministic system 
(r] = 0) possesses a stationary state uo, that satisfies 
f(u ,d x ,{fi}) = 0. Above a set of critical values {fj. c }, 
the system exhibits a spatial instability such that u reads 
u(x, t) = uo+e x ^ t+lkx u, where A (fc) is the linear growth 
rate and k the wave number of the instability. Close to 
{^c}, for A (fc) < (i.e. uo is stable), the profile of A 
already displays a maximum for a non null wave number 
close to the critical one k c . Hence, below the bifurcation 
point, when noise is present, among all the excited spatial 
modes, the one associated with the maximum growth rate 
(fc ~ fc c ) rules the dynamics. The dynamical behavior of 
the system is then characterized by pattern precursors. 
Figure Q] illustrates this property. It shows the instanta- 
neous S (fc, t) and the averaged S (fc) =< S(k, t) > t mod- 
ulus of the Fourier transform of the field u, also called the 
structure factor, for the supercritical Swift-Hohenberg 
equation 0] below the deterministic threshold. We can 
remark that the maxima of the function, both instanta- 
neous and averaged, already give the incoming critical 
wave numbers fc c = ±1 (FigQJ. 

To capture the dynamics of the system of Eq. 
(U) in an unified description close to the instability 
threshold, we introduce the ansatz u(x,t) = u D + 
A(T = et,X = yfex) e ik ° x u k + A (T, X) e~ lk - x u k + h.o.t 
0|, where A is the amplitude of critical mode, e is the bi- 
furcation parameter which is proportional to — fJ. c \\, 
Uk is the eigenvector associated to the critical mode, 
{T, X} are the slow time and spatial variables and the 
high order terms (h.o.t) are a series in the amplitude A. 
Introducing the above ansatz in Eq. {1} leads to the am- 
plitude equation that satisfies the solvability condition 



8 T A = eA- \A\ 2 A + d xx A + ^(X,T), (2) 

where £ (X, T) is a white Gaussian noise with null 
mean value and correlation (£ (X, T) £ (X 1 , T')) = 
and (i(X,T)i{X',T')) = 6(T-T)S{X'- X). The 
£ (X, T) term is then a linear combination of the original 



noise source term ((x,t). In absence of noise (r) = 0), 
when the bifurcation parameter e is negative, the system 
is characterized by a stable equilibrium state of ampli- 
tude A = 0, and when e becomes positive, the system 
exhibits a degenerate pitchfork bifurcation whose ampli- 
tude evolves as \A\ ~ \fe (deterministic value in Fig. [2]). 

The next step is to consider the global amplitude 
a (T) = A(T, X)dX/L, where L is the size of the 

transverse domain under study. This quantity is rele- 
vant to describe the bifurcation as soon as the trans- 
verse domain is either uniform or with very slow spatial 
variations of the control parameters. Experimentally, it 
will be achieved by choosing a restricted part L of the 
transverse pattern within which the control parameters 
are slowly varying. This restriction excludes all specific 
cases with fast spatial variations (e.g. local structures), 
which require specific analyses associated with the con- 
sidered transverse profiles, but they are out of the scope 
of this paper. The global amplitude a (T) then satisfies 
the Langevin equation 

d T a = ea- \a\ 2 a + ^{T), (3) 

where ?(T) = /^/ 2 £(T, X)dX/L is a Gaussian 
white noise with zero expectation value and correlation 
(? (T) f (T')) = 5 (r- T) . The general way of obtaining 
a solution of such Langevin equation is by use of a Fokker- 
Planck equation which provides us with a deterministic 
equation satisfied by the time dependent probability den- 
sity P(a, a; T) of the amplitude a [12] that reads 

d T P^d a [-ea+\a\ 2 a+^da)p + c.c. (4) 

where c.c. means "complex conjugate". The associated 
stationary probability density of the modulus of a is 



p s {\a\,£,r}) = Q{e,T))\a\e ^ , (5) 
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Figure 2: Bifurcation diagrams of the amplitude of the critical 
mode obtained for different solutions of Eq. ([2]). Namely, 
the deterministic value A, the expectation value |a maa; |, the 
linear and non-linear mean values (\a\)u nea r and (\a\) no „u„ ear 
respectively (see text for value definitions). 
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Figure 3: Stationary probability distribution P s (\a\) for dif- 
ferent values of the bifurcation parameter e. From left to 
right: e = -1; 0; 0.25; 0.5; 1. 7? = 0.2. 
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where Q{e,rj) = 2v / 2e~^/er/c(— -^=)^fwr\. The sta- 
tionary probability density is shown in Fig. [3]for different 
values of the bifurcation parameter e. The probability 
density function is not symmetrical with respect to its 
maximum so that the most relevant quantity for charac- 
terizing P g (|o|,£, J]) is its maximum and not its mean 
value as usually calculated e.g. in experiments. The 
value of \a\ corresponding to the maximum of P s (\a\, e, rj) 
occurs at the expectation value |a ma:E | given by 

\a>max\ = y 2 W 

When the bifurcation parameter e is driven far from 
(|e| 3> 1), this expectation value converge to zero 
(y/—T]/2e) for negative values of e and to y/e (as for 
the deterministic case) for positive values. Plotting this 
quantity versus s gives the supercritical bifurcation shape 
in presence of noise (dashed line in Fig. [2]). Let's now 
discuss the relevance of |a maa; | for describing the noisy su- 
percritical bifurcation. If we neglect the nonlinear term 
in the Langevin equation J3]), one can perform a linear 
stability analysis that provides us with the linear mean 
value of the amplitude modulus (\a\) linear = y—T^nf^. 
Note that this value diverges at the bifurcation point and 
is only valid for e < (Fig. [2]). The corresponding non- 
linear mean value {\<A) nonlinear ls computed numerically 
from the probability density P s (|a|). All the linear mean, 
non-linear mean, expected and deterministic values of the 
critical mode amplitude are reported on Fig. [2] for com- 
parison. The interesting region is located in the vicinity 
of the bifurcation point where the behaviors of the differ- 
ent curves strongly differ. The linear mean value and the 
deterministic value do not correspond to a realistic physi- 
cal behavior since the amplitude never diverges at thresh- 
old and we are considering a noisy system respectively. 
Only the nonlinear mean and expected values can mimic 
the supercritical bifurcation in presence of noise. How- 
ever, as we have mentioned earlier, due to the asymmetry 
of P s (\a\, £, 77), the most relevant value for describing the 
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Figure 4: Influence of the noise level r\ on the shape of the un- 
perfect supercritical bifurcation. The continuous curves cor- 
respond to the analytical expectation value |a ma2 ,| of Eq. (J6}. 
The squares and losanges are obtained from the numerical 
simulations of the stochastic Swift-Hohenberg equation Eq. 
©• 

evolution of the amplitude versus the control parameter 
is the expectation value |amax| - As this value includes the 
noise level 77, we have performed numerical simulations 
including noise to study the influence of noise on the bi- 
furcation shape. Figure S] clearly depicts the continuous 
deformation and evolution of the bifurcation shape with 
the level of noise. It also shows the very good agree- 
ment between the numerical values of |a maa; | obtained 
from the numerical simulations of the stochastic super- 
critical Swift-Hohenberg equation [3| and its analytical 
value [Eq. JfjJ]. Thus, the expression of |a ma 2;| is the 
relevant quantity for describing the noisy supercritical 
bifurcations including the noise level and the bifurcation 
point location (e — 0) of systems satisfying the Eq. {1} . 
Regardless of the sign of e, the width of the station- 
ary probability distribution decreases as |e| increases far 
from (Fig. [3}. Therefore, close to the bifurcation point 
(e = 0) the value of P s (|a m0 x|) exhibits a minimum ver- 
sus e (e « 0.25 in the case of Fig. [3]) which is shifted 
from this bifurcation point. The width of the stationary 
probability density is then maximum corresponding to a 
dynamics characterized by large amplitude fluctuations. 
This shift between the minimum of P s (\a max \) and the 
bifurcation point (e = 0) reminds the bifurcation post- 
ponements as in [15(. After a straightforward calculation, 
the value of the bifurcation parameter (e m i n ) associated 
with the minimum of P s (\a max \), for a given value of 
noise intensity, is 

(77) w 0.55^7- 

It allow us to propose a criterion for the localization of the 
bifurcation point: for a given noise intensity (77) the bifur- 
cation point is localized at e m in — 0.55^, where s m in is 
the value of the control parameter for which P s (\a ma x |, e) 
is minimum. 

In order to completely validate the universal law of 
Eq. ([6]), we have applied our analysis to experiments re- 
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Figure 5: Experimental bifurcation diagram for the optical 
feedback system. The dots are the experimental measured 
values of I d the continuous curve its fitted analytical 

value (Eq. [6). The dashed line is the corresponding deter- 
ministic bifurcation deduced from the fit. The inset figure 
is a schematic sketch of the experimental set-up. LC liquid 
crystal layer; M feedback mirror; F input optical field; B 
backward optical field; d feedback length. 



alized on a noisy ID transverse system known to exhibit 
a supercritical bifurcation at the onset of roll pattern 
formation. The system is a nematic liquid crystal slice 
subjected to optical feedback (cf. Fig, [p based on the 
well known feedback optical system [la Il7j|. The corre- 
sponding stochastic model reads 0| 

d t u= (d xx -\)u+\F\ 2 + R\e l ° d " (e txa F)\ 2 + ^}( (7) 

where u(x,t) stands for the refractive index of the non- 
linear nematic liquid crystal (LC) layer; t and x are the 
time and transverse space variables scaled with respect to 
the relaxation time r and the diffusion length I a; Ris the 
mirror intensity reflectivity, a = d/ko where d is the slice- 
mirror distance and kg is the optical wave number of the 
field. F is the forward input optical field, its transverse 
profile is accounted for using F(x) — Fq exp(— x 2 /w 2 ) 
for a Gaussian pump beam of radius w. C and rj are the 
noise source and level respectively as defined in Eq. {If. 
The Kerr effect is parametrized by x which is positive 
(negative) for a focusing (defocusing) medium. 

Eq. ([7]) is similar to Eq. (U) and the spatial variations 
of the pumping beam around its maximum are slow (less 
than 5% for a domain width L including 4 to 5 rolls) due 
to the high transverse aspect ratio (2wk c /2-K > 30), so 
that the previous analysis can be applied to describe the 
supercritical bifurcation of our experimental noisy sys- 
tem. On Fig. [5]we have plotted the experimental record- 
ings of the amplitude expectation value |a maa: | together 
with its analytical expression [Eq. |(6])] for a transverse 
domain width L including 4 to 5 rolls around the cen- 
ter of the Gaussian pumping beam. We can see that the 
analytical expression fits very well the experimental val- 
ues. It provides us with the deterministic threshold and 
noise level, I t h — \F t h\ 2 = 152 W.cm~ 2 and rj = 0.01 
respectively, rj is in very good agreement with previous 



determinations 13, 13]. So, the universal amplitude ex- 
pression of Eq. (J6j) is valid and relevant to describe the 
supercritical spatial bifurcation shape of our noisy sys- 
tem even and more generally to describe the supercritical 
bifurcations of ID systems in presence of noise . 

In conclusion, we have given an universal amplitude 
equation for ID systems in presence of noise. From this 
equation, we have derived an expression for the ampli- 
tude expectation value that fully describes the noisy su- 
percritical bifurcation. The agreement with experiments 
carried out for an ID pattern forming system is excellent. 
This amplitude equation can be applied to any second or- 
der transition of noisy temporal systems. 
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